Generalized survival in equilibrium step fluctuations 
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We investigate the dynamics of a generalized survival probability S(t, R) denned with respect 
to an arbitrary reference level R (rather than the average) in equilibrium step fluctuations. The 
exponential decay at large time scales of the generalized survival probability is numerically analyzed. 
S(t, R) is shown to exhibit simple scaling behavior as a function of system-size L, sampling time 
St, and the reference level R. The generalized survival time scale, r s (R), associated with S(t,R) is 
shown to decay exponentially as a function of R. 
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Controlling the stability of nanostructures is an im- 
portant fundamental issue in nanoscience. A key prob- 
lem in this respect is the random stochastic interface dy- 
namics associated with equilibrium nanometer scale ther- 
mal fluctuations unavoidabh/present in all nanosystems. 
Very interesting questions have recently been 

addressed regarding the complex dynamics of fluctuating 
interfaces in the context of first-passage statistics, which 
seems to be the appropriate framework in investigating 
the time it takes for a nanodevice to first fluctuate into an 
unstable state. It turns out that useful dynamical quanti- 
ties such as persistence probability P(t) 0,0 El ( Le - the 
probability that a stochastic variable does not return to 
its initial value over time t) and survival probability S(t) 
0,11] (i-e. the probability that a stochastic variable does 
not cross its average level up to time t) can be numeri- 
cally and experimentally investigated for interfaces with 
dynamics governed by various kinetic mechanisms (such 
as high-temperature attachment /detachment of atoms at 
the step edge or low-temperature step edge diffusion of 
atoms) to gain insight into the stability issue. 

Much work has been devoted 0,0 over the last decade 
in understanding equilibrium fluctuations on vicinal sur- 
faces, mostly using the dynamic scaling approach. If 
h(x,t) is the dynamical height (with respect to the cho- 
sen reference position which is defined to be the average 
position, i.e. the h = line) fluctuation of a thermally 
fluctuating step as a function of the lateral position x 
and time t (where t is also measured from an arbitrary 
time origin), then h(t) at each value of x is a stochas- 
tic dynamical variable by virtue of equilibrium thermal 
fluctuations. Because of the spatially extended nature of 
the step fluctuations through its dependence on both x 
and t, the problem is non-Mar kovian, and persistence Q 
and survival concepts should be particularly relevant 
statistical tools in understanding the complex problem of 
surface fluctuations. 

In this paper we introduce the new concept of a gen- 
eralized survival probability which enables us to probe 
deeper into the nature of the stochastic process of inter- 



face step fluctuations. The generalized survival proba- 
bility is the probability S(t, R) that a given lateral step 
position x with a height (i.e., step fluctuation measured 
from the equilibrium step position) h(x, t) at time t does 
never cross a pre-assisgned reference level of the height, 
R, throughout the entire evolution. The particular case 
with R — (i.e., the probability of the dynamical step 
height not returning in time t to its average ("equi- 
librium") R = level) has been studied recently || 
both analytically and experimentally, and it has been 
shown to exhibit an exponential decay at large times, 
S(t) oc exp(— t/r s ), where t s is the survival time scale 
that provides information about the underlying kinetics. 
The resulting surface step fluctuation survival probabil- 
ity S(t) = S(t, R = 0) and the associated time scale 
t s have also recently been studied experimentally using 
dynamical scanning tunneling microscopy (STM) on dif- 
ferent metallic systems: Al steps on Si (111) surface at 
high temperatures, and As and Pb (111) surfaces at rel- 
atively low temperatures ,9|. In this paper we show nu- 
merically that S(t, R) also has an exponential behavior at 
large time, S(t,R) cx exp(— t/r s (R)), where t s (R) is the 
generalized survival time scale. Our study reveals the de- 
pendence of t s (R) on the system size L, sampling interval 
8t, and reference level position R, allowing us to establish 
the complete scaling form of S(t,R). In particular, the 
sampling interval (i.e. the time between successive mea- 
surements) lOj turns out to be an essential ingredient in- 
herent in any real experimental measurement procedure. 
Also the study of the dependence of the generalized sur- 
vival time scale on the choice of the reference level R, 
which turns out to be exponential, should have partic- 
ular importance for understanding the effect of thermal 
fluctuations at the nanoscale. 

In this study we consider the case of the high- 
temperature step fluctuations dominated by atomistic 
attachment and detachment (AD), where the step edge 
is known to be well described by the coarse-grained 
second-order non-conserved linear Langevin equation, 
also known as the Edwards-Wilkinson (EW) equation 
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dh(x, t) 

at 



= V 2 h(x,t) + rj(x,t), 



(1) 



where V 2 refers to the spatial derivative (with respect 
to x, the lateral position along the step), r](x,t) with 
{rj(x,t)r](x , t )) = 2DS(x — x )8(t — t ) is the usual un- 
corrected random gaussian noise corresponding to the 
non-conserved white noise associated with the random 
AD process, and D is the noise strength. The AD pro- 
cess, thought to be extremely important for relatively 
high-temperature step fluctuations, has been extensively 
studied in the literature using the EW equation 0. 

For equilibrium step fluctuations, we define the gen- 
eralized survival probability with respect to the height 
reference level R, S(t,R), as the probability for the 
height variable to remain consistently above a certain 
pre-assigned value "i?" over time t: 

S(t, R) = Prob { h(x, t') > R, V t < t' < t + t }, (2) 

where h(x,t) is the dynamical height of the interface at 
a fixed lateral position x at time t, and to is the initial 
time of the measurement. Although the above definition 
involves the dynamical variable h(x,t) defined for a par- 
ticular lateral position x, we take a statistical ensemble 
average over all lateral positions to obtain a purely time 
dependent stochastic dynamical quantity S(t,R). Ob- 
viously, another quantity that can be measured is the 
probability for the height stochastic variable to remain 
below the reference level up to time t. Since in our case 
the dynamics of the interface fluctuations obeys a linear 
stochastic equation, the interface preserves the up-down 
symmetry along the direction perpendicular to the step 
edge. As a consequence, in what follows we consider the 
average of the probabilities of remaining always above R 
and below —R, with R > 0. 

The generalized survival probability function, S(t,R), 
defined in Eq. J2J above, leads to a hierarchy of gen- 
eralized survival time scales, t s (R), if the steady-state 
decay of S(t, R) in time follows an exponential trend, 
S(t,R) ~ e"*/ 1 " 3 ^. As we show below, this indeed is 
obtained for Edwards-Wilkinson equilibrium step fluc- 
tuation phenomena, allowing us to define and measure 
the non-trivial survival time scale t s (R), < R < R ma x, 
that varies between t s (R — 0) and T s (R max ), where 
t s (R = 0) is the usual survival time scale and T s (R max ) 
is the survival time with respect to the highest refer- 
ence level R m ax that can be defined for a model with 
finite roughness (i.e. rms fluctuations of the height vari- 
able with respect to the average) . R m ax is limited by the 
maximum value of the height fluctuation amplitude. Ob- 
viously, S(t,R) and t s (R) are natural generalizations of 
the survival probability S(t) and the survival time scale 
t s , respectively, to the more complex concept of distribu- 
tion of generalized survival times with limiting behavior 
(i.e. R — 0) providing the usual survival time. 



The exponential decay at large time of S(t, R) that 
we find numerically is not surprinsing. The generalized 
survival probability with respect to the reference level R 
can be regarded as the probability Z(t) of no zero cross- 
ing of the new stochastic variable H(x,t) — h(x, t) — R. 
What we are looking for is the probability for the stochas- 
tic variable H(x,t) to remain positive up to time t 
(or, equivalently, the probability for h(x,t) + R to re- 
main negative over time t). This type of question for 
the gaussian stationary processes with zero mean has 
been addressed by mathematicians for a long time [l^ . 
The no zero crossing probability is traditionally investi- 
gated in conjunction with the autocorrelation function, 
Ch = {H(x,tx)H{x,t2)) (where (...) represents an aver- 
age over all realizations of H(x.t) arising from the ther- 
mal noise source). It is known jl3j that for a stationary 
gaussian process (i.e. Ch = f(\ti — tz\) = f{t)) with 
an autocorrelation function decaying faster than \jt at 
large t, the asymptotic behavior of the no zero crossing 
probability is exponential, Z(t) oc exp(— fj,t). The auto- 
correlation function Cn(t) itself has been shown [5j to 
be stationary at late times and to decay exponentially. 
This, along with the exponential decay of Z(t), ensures 
an exponential decay for S(t, R). 

In order to numerically simulate the process described 
by Eq. (JJJ, we have used discrete stochastic Monte 
Carlo simulations of the corresponding atomistic solid- 
on-solid model, the extensively studied Family model [|J, 
which belongs asymptotically to the Edwards- Wilkinson 
universality class [ll|. The Family model in (1+1) 
dimensions (i.e. one spatial variable and one temporal 
variable) is characterized by (3 = 1/4, a = 1/2 and 
z = a//3 = 2 |(|, where the growth exponent (3 is the 
rate of change of interface width (or roughness) in the 
transient regime (w(t) ~ t@), the roughness exponent a 
shows the saturation of the width for a system with fixed 
size L in the steady state regime (w(L) ~ L a ) and z is the 
dynamical exponent. This model involves the traditional 
random deposition (at a rate of one complete monolayer 
during one unit of time) and surface relaxation such that 
the adatoms are searching for the sites with the minimum 
local height. We have taken the relaxation length to be 
the lattice constant and we have applied the usual peri- 
odic boundary conditions. Typical sizes (i.e., number of 
lattice sites) used in this numerical work are 100 — 900, 
and the averaging procedure implies a number of at least 
10 5 independent runs. All the measurements correspond 
to the steady state regime where the interface roughness 
has reached a time independent equilibrium value (i.e., 
to ^> L z in Eq. (0)). We also mention that the small- 
est value for the sampling time is 1. We emphasize that 
our use of Family model is just a matter of convenience in 
simulating the EW equation [£j ; our results are simply an 
exact discrete stochastic simulation of the EW equation. 

Our results for the generalized survival probability and 
the associated time scale are presented in Figs. ^ and 
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FIG. 1: The generalized survival probability, S(t,R), for the 
discrete Family model. The dashed lines are fits of the long- 
time data to an exponential form. The system size is L = 100, 
the sampling time is 5t = 1.0 and the reference level R takes 
four different values: 0, 1, 2 and 3 (from top to bottom). The 
inset shows the dependence of the generalized survival time 
scale Ts(R) on the reference level value (up to R = 5). The 
continuous curve represents a fit to an exponential decay of 
the generalized time scale vs. R. 



121 S(t, R) is simply computed as the fraction of sites 
which, starting above (below) the level R (-R) at time 
to, have not crossed the reference level up to a later time 
to + t. In Fig. ^ we show that, as expected, the gen- 
eralized survival with respect to an arbitrary reference 
level R follows an exponential decay at large times. The 
only varying parameter in Fig. ^ is the reference level 
R. We have considered six values for R, R = 0,1,..., 5 
(only the first four curves are displayed due to the limi- 
tations imposed by the quality of the statistics, since as 
R increases it is less probable to have a reasonable num- 
ber of lattice sites with height variables above (below) R 
(— R)). The dashed lines are fits of the long-time data 
to an exponential form, S(t,R) oc exp(— t/r s (R)). The 
upper curve has R = and corresponds to the usual 
survival probability previously studied in Ref. 5]. How- 
ever, all the other curves are new and they prove that the 
generalized survival probability decays exponentially in 
the long-time limit, with an associated time scale, t s (R), 
which decreases with the reference level value. As shown 
in the inset of Fig. ^ the dependence of r s (R) on R is 
exponential, but clearly more work is needed in order to 
understand this trend. 

In Fig. |2 we have used several lattice sizes, sampling 
times and reference levels in order to identify the scaling 
behavior of S(t,R). In panel (a) we show the general- 
ized survival with respect to level R=l, measured using 
5t = 1, for two system sizes: L = 100 and L = 200. We 
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FIG. 2: The generalized survival probability, S(t,R), for the 
Family model. The dashed lines are fits of the long-time data 
to an exponential form. Panel (a): L = 100 (lower curve) 
and L = 200 (upper curve), using fixed sampling time St = 1 
and reference level R = 1. Panel (b): St — 1 (lower plot) 
and St = 16 (upper plot) with fixed system size L — 200 
and reference level position R = 1. Panel (c): Scaling of 
S(t, R) using three different system sizes with sampling times 
and reference levels varied such that St/L z and R/L a are 
kept constant (i.e. St/L z = 1/10000 and R/L a = 1/VW0). 
A perfect collapse of the curves with L=100, 400 and 900, 
respectively, occurs when using z = 2.03. 
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observe that the underlying survival time scale increases 
rapidly with L. In fact, t s (R) for a fixed R is expected to 
grow proportionally to L z |5|. However, we obtain that 
t s (R = 1) ~ 103 for L = 100, and t s (R = 1) ~ 429 
for L = 200, so the measured generalized survival time 
exhibits a small deviation from the expected value of 
103 x 4 = 412. We find that this small effect is due to the 
dependence of the generalized survival on sampling time 
St. This is clearly seen in panel (b). It turns out that 
a system with a fixed size (L = 200) is characterized by 
different values of r s (R) if the sampling time of the mea- 
surement is adjusted. We observe that t s (R) increases 
weakly as the sampling time is increased. One might ar- 
gue that this effect is very small and could be neglected, 
but we have found that the effect of the sampling time 
on the measured generalized survival probability has to 
be taken into account in order to find the complete scal- 
ing function of S(t,R). In addition, this effect is even 
stronger for systems with slower dynamics (i.e. larger z) 
Interestingly enough, we note that fixing the refer- 
ence height level in the generalized survival probability 
problem introduces an additional length scale, that is re- 
lated to the steady state value of the interface width, 
i.e. L a . Indeed, in panel (c) we look at three different 
systems with L — 200, 400 and 900, respectively, and 
the generalized survival curves are calculated for R = 1, 
2 and 3, respectively, i.e the level R is varied propor- 
tionally to L a , with a = 1/2 as appropriate for the EW 
equation. In addition, the sampling time for each of these 
three cases is also varied, St <x L z (z = 2), so we have 
considered St = 1 for L = 100, St = 16 for L = 400, 
and St = 81 for L = 900, respectively. A perfect col- 
lapse of the curves S(t,R) vs. t/L z occurs when using 
z = 2.03, which agrees with the expected value z = 2, 
characteristic for the EW dynamics. 

This numerical analysis allows us to conclude that the 
scaling form of the generalized survival probability is 

S(t, L, R, St) = f{t/L z , R/L a , St/L z ), (3) 

where the function fix, y, z) decays exponentially for 
large values of x. The rate of this decay decreases rather 
rapidly as y is decreased and increases rather slowly as z 
is decreased. Note that for y = we recover the scaling 
form of the usual survival probability with R = || . 

To conclude, we have shown that the generalized sur- 
vival probability of equilibrium step fluctuations on vic- 
inal surfaces with Edwards-Wilkinson dynamics decays 
exponentially at long times. We have investigated the 
associated generalized survival time scale that depends 
on the system size L, sampling time St, and the choice 
of the reference level R. In particular, the depedence 
of t s (R) on R, which based on our preliminary investi- 



gations seems to have an exponential trend, should be 
useful in understanding the stability of thermally fluctu- 
ating interfaces. We have also shown that the generalized 
survival probability exhibits simple scaling as a function 
of L, St, and R. Our numerical results on S(t, R) can be 
easily extended to fluctuating interfaces characterized by 
different dynamical evolutions (such as low-temperature 
step edge diffusion limited kinetics) belonging to differ- 
ent universality classes. Our goal here, using the example 
of the step fluctuations process characterized by the EW 
universality class, is to establish the generalized survival 
probability as an important statistical concept in study- 
ing thermally fluctuating interfaces. 

Finally, we mention that the generalized survival prob- 
ability could be experimentally measured using dynam- 
ical STM step fluctuations data, opening the possibil- 
ity for a direct approach to the crucial issue of interfa- 
cial stability. Our theoretical considerations for S(t, R) 
should also be useful in understanding the dynamical 
evolution of other physical processes [14| where a first- 
passage statistics has proven to be an useful concept. 
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